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1 Introduction 



In the previous chapters of this book, experimental results and the theory of elastic and inelastic 
neutron scattering have been presented under the assumption that the magnetic moments of the 
neutrons are randomly oriented. It has been shown that details about the physical properties of 
a system are extracted by analysing the momentum and the energy of the scattered neutrons. It 
is intuitively imaginable, however, that measuring the spin state of the neutron after scattering 
relative to its state before the scattering process should provide us with additional information. 
To that end the cross-sections for neutron scattering must now also take into account the 
relationship between the spin of the neutron with the physical properties of the target. 

The polarization of a neutron is defined as 

P = 2<3) = <*) > (1) 

where cr are the Pauli matrices. Clearly, |P| is equal to for a completely unpolarised beam 
and |P| = 1 if the beam is totally polarised. For intermediate values, the neutron beam is not 
in a well defined state and the spin part of the neutron wave-function must be described by a 
more general form \ — U X] + v Xi with \u\ 2 + \v\ 2 = 1. That is \u\ 2 and \v \ 2 are the probabilities 
that the neutron spin will be up or down, respectively. If a matrix operator p is defined like 

P = XX*=[ { J H2 J = 2(I + P-a), (2) 

then the polarization of the neutron is described by a three-dimensional vector wit components 
P = (23l(u^v ), 25s (vlv ), \u\ 2 — \v\ 2 ) ]IJ. The polarization of a neutron beam is accordingly 
defined as P = ■h ^ Pj where N is the total number of neutrons and the sum runs over the 
polarization vector of the individual neutrons j. 

The cross-section a and the polarization of the scattered beam Pf can be expressed as a 
function of the density matrix p, the polarization vector of the incident neutrons Pj and the 
interaction potential v between the target and the neutron. In its most general form, Pf — 
Trpv^&v/Trpv^v. 

Neutron scattering with polarized neutrons has been used in fundamental and condensed matter 
physics since many years despite the low flux of polarized beams. One of the first applications 
was the study of spin density distributions in ferromagnets, following the pioneering work of 
Shull and Nathans . In those days the polarization of the scattered neutrons was established 
by measuring their transmission through a magnetized block of iron. As predicted by Halpern 
and Johnson |3j the polarization of the scattered neutrons, Pf, depends on the orientation of 
the scattering vector Q with respect to the polarization of the incident neutrons, Pj, like 

P f = -Q(Q ■ Pi) (3) 

where Q = Q/|Q|. In other words for Pj || Q all magnetic scattering is spin flip. Therefore 
polarization analysis in neutron scattering provides an excellent method to distinguish between 
nuclear and magnetic scattering. In 1969 the classic paper by Moon et al. [|J appeared that 
explains in simple terms one- dimensional polarization analysis (nowadays called longitudinal 
polarimetry) of neutrons for elastic as well as inelastic neutron scattering. They demonstrated 
the polarization dependence of the nuclear and magnetic scattering. 
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Another, rather different application of polarized neutrons is their use for attaining extremely 
high energy resolution by measuring changes in the neutron beam polarization caused by in- 
elastic scattering. In 1972, Mezei || suggested to use the precession of the magnetic moment of 
the neutrons in a magnetic field as an internal clock. By means of the so-called neutron spin- 
echo technique energy resolutions of the order of nano-electron volts can be achieved enabling 
the investigation of slow dynamics, for example in the critical region of magnetic systems or in 
polymers and glasses. 

Nowadays, polarized neutron scattering is a fast developing experimental method that finds 
applications in various fields of condensed-matter research. Examples are 

• determination of magnetic structures and spin densities, 

• identification of magnetic fluctuations and their different modes, 

• separation of coherent from incoherent processes. 

In the following, necessarily incomplete sections, we shall provide a presentation of the polariza- 
tion dependence of neutron cross-sections and show how the different scattering processes can 
influence the polarization of the neutron beams. We will then explain how polarized neutron 
beams can be produced and the polarization determined after scattering. Finally we shall give 
examples, where the technique of polarized neutron scattering can provide new insight into 
physical processes in condensed matter research. For a more detailed introduction into the field 
of polarized neutron scattering with refer the interested reader to the book of Williams |J. 

2 Elastic neutron scattering cross-section for polarized 
neutrons 

The theory of elastic neutron scattering taking into account polarization effects has been derived 
by Blume in 1963 0. The complete description of the scattering process involving both nuclear 
and magnetic interactions can be given by means of two master equations. The first one gives 
the total neutron cross-section which depends on the polarization P« of the incident neutron 
beam as follows: 

a = NN* + D± ■ Dl + Pi(D±N* + B* ± N) + iP f (Dl x D ± ), (4) 

where a is the total cross section expressed in [cm 2 ]. For simplicity, the contribution of the 
nuclear spins is neglected. iV = iV(Q) = ^ h exp (iQ ■ r») is the structure factor of the atomic 
structure that depends on the scattering vector Q and the scattering lengths of the individual 
nuclei bi\ Dj_ is the magnetic interaction vector with = Dj_(Q) = Qx (p(Q) X Q)- p(Q) 
is the Fourier transform of the magnetic moment distribution and Q = Q/|Q|. Therefore, 
only magnetic components perpendicular to the scattering vector participate in the scattering 
process. The scalar of the polarization vector P; reflects the degree of polarization of the 
neutrons, being equal to ±1 for a fully polarized beam. 

Eq. H shows that the neutron cross-section depends only on the square of the chemical and 
magnetic structure factor if a non-polarized neutron beam (Pj = 0) is used. For a fully 



A 



polarized beam (|P,,| = 1), two additional terms contribute to the scattering, namely the 
magnetic-nuclear interference term and the chiral term, respectively. The magnetic-nuclear 
interference term being proportional to Dj_iV* + D^iV yields only a non- vanishing contribution 
to the neutron cross-section if a Bragg reflection is due to nuclear and magnetic scattering, like 
in ferromagnets and in non-centrosymmetric antiferromagnets with propagation vector Qo = 0. 
The chiral term x is non-zero whenever is not parallel to D^_, as it is the case e.g. 
for a helicoidal magnetic structure. 

The second master equation provides the polarization of the neutron beam after the scattering 
process relative to the polarization of the incident neutron beam: 

P/tr = PiNN* 

+ (-l)Pi(D ± • Dl) + D ± (P, • Dl) + Dl(Pi • D ± ) 

+ B ± N* + ~DIN + i(D ± N* -~DIN) xPi 

+ «D ± x Dl, (5) 

where P/ is the polarization vector of the scattered neutrons. On the one hand, Eq. |5] shows 
that pure nuclear scattering (Dj_ = 0) leaves the polarization of the neutron beam unchanged. 
On the other hand, polarisation of the scattered beam is obtained either, as we will see in the 
next chapter, by scattering neutrons on mixed nuclear-magnetic Bragg reflections or from a 
helicoidal magnetic structure. In the latter case, with Dj_N* + D^N = 0, a polarized beam 
with a polarization given by 

P/ = — ; — = D ± D1 (6) 

is created. We point out that a measurement of the chiral term provides the helicity of a 
helicoidal magnetic structure as has been shown by Shirane et al. ||. In the general case, the 
polarization vector of the neutron beam after scattering is rotated with respect to Pj and its 
length is not necessarily equal to |Pj|. The term "polarization analysis" therefore refers to the 
determination of the direction and length of P/. 



3 Production of polarized neutrons 



For a measurement of the polarization dependence of cross sections various techniques to pro- 
duce and analyze polarized neutron beams have been developed. Depending on the required 
phase space properties of the beams, i.e. continuous vs. pulsed, energy, divergence, type of 
detector, etc., different methods for the spin analysis are used. The most common methods 
are diffraction from single- crystal polarizers (mostly Heusler), reflection from magnetized thin 
film multi-layers or supermirrors, and absorption of the non-wanted spin state by means of 
polarized 3 He. A recent review of these techniques can be found in Ref. ||. 



3.1 Single-crystal polarizers 



This method produces a polarised neutron beam by taking advantage of the magnetic-nuclear 
interference term in Eq. f| and Eq. [j| If a magnetic field is applied to a centro-symmetric crystal 
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so that all the magnetic moments are saturated and aligned perpendicular to the scattering 
vector Q, the neutron scattering cross-section for Bragg scattering is given by (set P; = and 
D ± = in Eq. |) 

a = N 2 + D\. (7) 
The second master equation, Eq. [|, yields for Pj = 

2iVD ± 2N~D ± 



f 



N 2 + D 2 ' 



(8) 



Hence, the diffracted beam from a single-crystal is completely polarized if there is a Bragg 
reflection with |Dj_| = \N\. Typical examples are the (111) reflection of Heusler CU2M11AI (d- 
spacing 3.43 A) and the (200) reflection of the alloy Coo.92Feo.08 (<i-spacing 1.76 A). Other single 
crystals like Fe 3 4 or Fe 3 Si have also been considered but are less used. All these crystals can 
be used to produce polarized and monochromatic neutron beams and to analyze the energy and 
polarization of neutron beams. Therefore, single-crystal polarizers are used for single- crystal 
diffractometers and in triple-axis spectroscopy. Depending on the requirements on neutron 
energy and resolution, different <i-spacings must be considered. Recently, the quality and the 
reflectivity of Heusler monochromators has been improved considerably |L0| that will allow to 
use these crystals at relatively short neutron wavelengths. 



3.2 Thin films 



Total reflection from magnetized thin films can be used to produce polarized neutrons. The 
angle of total reflection for a ferromagnetic film is given by 

9f = Xy/N(b±p)/ir, (9) 

where A is the neutron wavelength. N is the nuclear density, and b and p the nuclear and mag- 
netic scattering lengths, respectively. Thus, by an appropriate choice of materials, a polarized 
beam can be produced by total reflection. For the special case b = p all reflected neutrons are 
polarized. Unfortunately, the reflection angles are only reasonably large for cold neutrons: For 
example, Fe 5 oCo 48 V 2 has b ~ p and 9 C ~ 0.4° for A ~ 4 A fTT| . 



The angles of reflection can be significantly improved by adding artificial magnetic and non- 
magnetic layers that reflect neutrons at small angles above 8 C . Such artificial multi-layers 
(supermirrors) have been produced by physical vapor deposition by Mezei for the first time 
fl2f . Typical materials combinations are Co/Ti, Fe/Si, and Fe 5 oCo4 8 V 2 /TiNi x 0, [13], [14], [lj| [L6 



The latter combination exhibits a remanent magnetization and can therefore be used as a spin 
selective device, where no spin flipper is necessary anymore | I7|| . Recently, the number of layers 
has been increased steadily thus leading to reflection angles for polarized neutrons of the order 
of 0.3° for A ~ 1 A. These modern devices can now also be used as white beam polarizers for 
thermal neutrons. 



3.3 Spin niters 



A major drawback of polarizing single-crystals and thin films is their decreasing efficiency with 
increasing neutron energy, i.e. short wavelength, and the small divergence that they accept (see 



Table |T]). On the other hand they are maintenance free and easy to use. Therefore, polarizing 
filters with broad-band characteristics and minor restrictions on divergence are of significant 
interest for neutron scattering in particular for pulsed spallation sources and spectrometers 
with large area-detectors. 

Table 1: Performance and applications of various neutron polarizers. The quoted values are only 
approximate. The notation is TAS: triple-axis spectrometer, DAX: double axis diffractometer, 
NSE: neutron-spin-echo, REF: reflectometer, TOF: time-of-flight, SANS: small angle neutron 
scattering. The quality factor is defined by Q = TP 2 , where the transmission/reflection T and 
polarization P are taken from the literature. 



technique 



beam .E-range (meV) 



instruments 



T 



P Q ref. 



Heusler 
super- mirror 
3 He 
p targets 
SmCo 5 



fixed A 
white 
white 
white 
white 



E < 80 
E < 20 
E < 2000 
E > 2000 
20 < E < 180 



TAS, DAX 0.62 0.95 0.56 

NSE, TAS, REF, TOF 0.9 0.95 0.81 

DAX, TAS, TOF, SANS 0.4 0.80 0.26 

SANS 0.4 0.80 0.26 

not implemented 0.3 0.75 0.17 



El 




Hi 




Whereas Heusler and supermirrors are well established but still progressing techniques, the 3 He 
spin filters using direct optical pumping of meta-stable 3 He have improved during the last 
few years [ TE\ . The basic idea behind the filter-technique is the polarization dependence of the 
transmission that can be written in its most simple form as || 



T(A) = exp (— a Nd) cosh(a p Nd), 



(10) 



where A is the neutron wavelength, d the thickness of the filter and N the 3 He density. <To and 
cr p are the spin-independent and polarization dependent cross-sections of 3 He, respectively. If 
the filter is not perfectly polarized, a significant portion of the correctly polarized neutrons will 
be absorbed. Because the absorption increases with increasing A, the thickness d of the filter 
must be optimized for the wavelength band to be used even if o"o is small. 

3 He filters are now in regular use at the ILL on several instruments. Due to wall relaxation 
of polarized He nuclei, the polarization and transmission of the filters decreases with time and 
they have to be exchanged almost daily. It is foreseen that 3 He filters will further improve 
and find applications in particular at pulsed neutron sources and instruments with large area 
detectors. The method of using spin exchange of 3 He with optically pumped Rb vapor 



is 



progressing too and may challenge the meta-stable type of pumping. One major advantage of 
the latter technique is that the filter has not to be exchanged during an experiment. 



The development of SmCos polarizing filters has been conducted at ISIS. If the problems of 
depolarization of the neutrons within the filter and of 7-heating in intense neutron beams can 
be solved one may obtain a quality factor of Q ~ 0.25 ||T^| that is lower than the maximum to 
be achieved for a future 3 He filter. 



Polarized hydrogen can also be used as polarizing filter. Moreover, the spin dependent interac- 
tion of the neutrons with the protons can be used for contrast variation of hydrogen-containing 
materials in small angle neutron scattering experiments [22]. We defer the interested reader to 
the literature. 
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4 Determination of Form Factors and Spin Densities 



Polarized neutrons allow to measure magnetic densities with improved accuracy as compared 
to standard diffraction methods. The method presented below applies to magnetic structures 
described by a propagation vector Q = 0. For a paramagnet, a ferromagnetic component 
can be induced by applying an external magnetic field. According to Eq. ^ in the case of 
mixed nuclear-magnetic Bragg reflections with real structure factors, the intensity ratio R of 
scattered neutrons polarized by an external magnetic field along the +z or — z direction, where 
z is perpendicular to the scattering plane (that contains Q), is given by 

I +z _ N 2 + 2ND Z ± + D z2 _ (N + D z ± ) 2 
~ T~ z ~ N 2 - 2ND Z ± + Df ~ (N - Df) 2 ' ( ' 

D\ is the projection of the magnetic interaction vector Di along the z-axis. Determination of a 
spin density with polarized neutrons consists of measuring R at many different Bragg reflections 
(hkl). As the crystal structure and hence the chemical structure factors N are presumably 
known (N depends on the Miller indices), the method provides usually directly the values 
of the magnetic structure factors. For small magnetic amplitudes, polarized neutrons give an 
enhanced sensitivity compared to unpolarized neutrons that yield an intensity I np = N 2 + \D\- 
The factor 2/3 comes from the spherical averaging of Di in the second term of Eq. |] over all 
directions with respect to Q for a non-magnetized isotropic sample. 

Namely, considering a typical example with D± = 0.1N yields I np ~ l.OliV 2 , while the contrast 
as measured with polarized neutrons, R = 1.21N 2 /0.81N 2 = 1.49, is rather large. There- 
fore polarized neutrons are particularly well suited for measuring maps in compounds with 
small magnetic moments (example: heavy fermion systems) or with seriously diluted magnetic 
moments (example: molecular magnetic crystals). 



To extract the magnetic moment density from the data, a classical Fourier calculation is usually 
performed p3| . Because 

P(Q) = / / / m(r) exp(zQ ■ r)rf 3 r, (12) 



one can obtain the spin density m(r) in real space by the inverse Fourier transform through 
the relation 

™(r) = -^5>(Q)ex P HQ-r). (is) 

Q 

As the cloud of unpaired electrons that are responsible for magnetism is extended in real space, 
the magnetic form factor decreases with increasing Q and equivalently with increasing Bragg 
indices (hkl). To obtain precise measurements, data are to be taken up to large values of scat- 
tering vectors Q. To that end short-wavelength neutrons are to be preferred, and instruments 
dedicated to such measurements provide usually hot neutrons, like the 2-axis diffractometers 
D3 at the ILL and 5C1 at the LLB. The layout of such an instrument is presented in Fig. [I| 
Until now, no polarized neutron diffractometers have been built at a spallation source. In any 
case, however, the data set is restricted to finite values of h, k, I, which leads to oscillatory 
distortions in the spin density maps m(r) due to finite size effects ||23|| . Therefore, other re- 
construction methods of the spin density map have been developed. The two most often used 
techniques are either based on information theory like the maximum entropy method p|| or 



on the multi-polar expansion of the electronic density |25[ . The latter method models the spin 



density by a set of parameters that have to be determined by standard least-square fitting 
calculations. 

Figure 1: Schematic arrangement of a two- axis spectrometer used for the determination of 
magnetic densities. The monochromator produces a polarized beam with the neutron polar- 
ization perpendicular to the scattering plane. A small guide field prevents the neutron beam 
to depolarize. A spin-flipper allows to reverse the neutron polarization by 180° and hence to 
measure the flipping ratio R. A magnetic field saturates the magnetic moments of the sample 
along the neutron polarization. 



4.1 Magnetic Form Factors 



Ever since the experiments of Shull and coworkers in the 1960's in 3d ferromagnets p6|| , the main 
motivation to measure spin densities has been to gain a better insight in electronic distributions 
in solid-state materials. Since that time, measurements have been extended to paramagnetic 
metals and to the 4/-electrons in rare-earth compounds p7| . pSfl . 



In 3d-ferromagnets whereas the atomic form factor can be very precisely reproduced from 
spin-polarized Hartree-Fock calculations p9| , there is a strong indication from form factor 
measurements of a negative spin polarisation between the atomic sites in both Fe and Co. 
Also, the magnetic moment density-map in Ni |3(| and Pd |yj shows an aspherical d-electron 
distribution plus an orbital contribution. In order to allow for these effects, the form factor is 
usually written as 

/(Q) = ^(l + a)[<jo> + l-~, 



9 
+ 



1 + a) 
• 9 



' f orb 



a8{Q) 



(14) 



where g is the Lande factor; a is a parameter describing the fraction of negative spin polar- 
ization; 7 is the percentage of electrons in E„ orbitals which takes into account the orbital 



contribution ^0|. Ahu is a geometrical factor; (j ) and (74) represent the spherical and as- 
pherical part of the form factor, respectively. A comparison of the calculated and measured 
form factors for Ni yielding a uniform negative contribution equal to — 0.0091yUe/A 3 is shown 
in Fig. [| The spin magnetic moment per Ni-atom is fi sp i n = 0.656/is-B and the orbital contri- 
bution florbital = 0.055/is. 



Figure 2: Magnetic moment distribution of Ni in the [100] plane (taken from Ref. [50] 



On the contrary to metallic compounds with ci-electrons, /-electrons are well localised around 
the nuclei which allows to perform atomic calculations to obtain the spin distribution. As the 
orbital moment is usually different from zero, there is an significant contribution of the orbitals 
to the magnetic density. Also, the spin-orbit coupling is important which results in a mixing of 
the atomic wave functions. The form factors for the atoms of the rare-earth and actinide series 
have been calculated by Desclaux and Freeman |32[ using the relativistic Dirac-Fock theory. It 
was shown that the atomic form factor can be expressed as 

/(Q) = O'o) + P2O2) + C4O4) + c 6 (j 6 ), (15) 

n 



where 



oo 



0'i(Q)> = / U\r) 3l (Qr)A^dr. (16) 



o 



U(r) is the radial wave-function for the unpaired electrons in the atom, and ji(Qr) the Bessel 
function of Zth-order. The coefficients c, are tabulated in e.g. Ref. f33j. 



Among the rare-earth elements, Samarium represents a particular case as the orbital and spin 
contributions to the magnetisation almost cancel out so that the magnetic density map contains 
both positive and negative regions. This leads to a form factor which has a maximum located 
at a position different from Q = 0. The Sm form factor measured in SmCos is shown in Fig. [| 
An interesting effect is, that as the first excited crystal-field states are located at relative low 
energies, they become populated when the temperature approaches 300 K. Consequently, the 
magnetic moment of Sm is only fi ~ 0.04/i^ at room temperature and increases to fi ~ 0.38yU# 



for T = 4.2 K [34|. For all temperatures, however, the form factor of Sm has a strong orbital 
character. 

Figure 3: a) Experimental form factor for Sm at T = 4.2 K. The line corresponds to a calculation 
including crystal field, exchange and spin-orbit effects, b) Same for T = 300 K (taken from 
Ref. M). 



Finally, we should point out that a study of the spatial distribution and temperature dependence 
of the spin density allows to probe the spin susceptibility x(Q,0). As such the method can 
been used to investigate the nature of the electrons e.g. in superconductors. For example, in 



VaSi |35| it was shown that the spin susceptibility of the V electrons disappears upon entering 
the superconducting phase which is an indication of spin-pairing. On the other hand, no 
similar effect could be observed in the new heavy- fermion superconductors UPt3, UBei3 and 
CeCu 2 Si 2 . For the latter compounds, the spin susceptibility is temperature independent in the 
superconducting phase fl3~6"| . These results are of particular importance as they impose severe 



restrictions on the possible pairing mechanisms that can give rise to the electron pairing in 
these unconventional heavy-fermion superconductors. 



4.2 Magnetisation distribution in molecular magnets 

Molecular magnetism is a fast growing field in material science with potential important tech- 
nological applications in electronic devices. By building blocks of molecules which contain 
magnetic centres, magnetic interactions can be tuned and the aim is to synthesize organic com- 
pounds which exhibit magnetic ordering at room temperature. To that end, as the number of 
combinations offered by organic chemistry is almost infinite, it is essential that the mechanism 



of magnetic couplings originating from 2p electrons to be well understood [37 



In contrast to ionic systems where the electrons which carry magnetism are well localised 
around the nuclei, the magnetic density of organic compounds is distributed over all molecules 
due to covalency effects. The effect of delocalisation is even more pronounced when there is no 
magnetic ion in the molecule and magnetism is due to 2p electrons only [Q . Polarised neutron 
diffraction yields directly the distribution of electrons responsible for magnetism in organic 
materials which in turn can be directly compared to theoretical calculations for the electronic 



i n 



wave-functions and the chemical bonds |p9| . Spin density in molecular compounds can also be 
used to trace exchange pathways through the molecules, like when spin polarisation is found 
on atoms which in principle are non-magnetic. This is for example the case in the free radical 
nitronyl nitroxides NitPy(C=C-H). NitPy(C=C-H) builds zig-zag chains linked by C=C-H- • -O 
pieces where the hydrogen bridge two molecules. Significant spin population is found at the 
hydrogen positions (// ~ 0.04/z#) which indicates that the hydrogen bond is involved in the 
ferromagnetic exchange interactions between the molecules £I(| . 



Figure 4: Spin density projection in NitPy(C=C-H). The contour step for the pyridine cycle is 
0.008/ib/A 2 whereas for the Nit cycle a step is equal to 0.04/zb/A 2 (taken from Ref. j40fl). 



A typical example of spin delocalisation is found in the compound MnCu(pba)(H 2 0)3 • 2H 2 0, 
with pba=l,3-propylenebis(oxamato). The Mn 2+ and Cu 2+ ions are connected by oxamato 
bridges and build ferrimagnetic chains. The Mn 2+ ions carry a spin Smh — § and the Cu 2+ 
have an effective spin Sc u — §■ Therefore we have the situation where two magnetic metals are 
linked by organic species. The magnetisation, as obtained from polarised neutron diffraction 
|f41| , |44j], shows a positive spin population (i.e. the induced magnetisation is aligned along 
the applied magnetic field) for the Mn spins whereas it is negative for the copper magnetic 
moments. This reveals the antiferromagnetic nature of the intra-chain coupling. Interestingly, 
an important contribution to the spin density map is found on the neighboring oxygen and 
nitrogen ions and on the two central carbon atoms. Summing up the positive and negative 
spin polarisations individually, one obtains 5.1/ib and — 1.0/zb, respectively, which shows that 
the metallic ions have distributed their spin densities on the molecule. The magnetic moments 
distribution for MnCu(pba)(H 2 0) 3 • 2H 2 is shown in Fig. ||. 

Figure 5: a) Experimental spin density map of MnCu(pba)(H 2 0)3 ■ 2H 2 0. The contour step is 
0.005/! B /A 2 . The cont inuous line represents the positive spin distribution while the dotted line 
describes the negative magnetisation, b) Calculated spin density map for an isolated molecule 
with the DMol 3 method (after Ref. |4j]). The lowest contour is at 0.005/i B /A 2 



Fig. ||also shows the theoretical spin density for the 'CuMn' molecule projected onto the oxamid 
mean plane. The theory is based on the local-spin density-functional principles of Perdew and 
Wang fl42f . Calculations for the cation (2+) in vacuum are done with the DMol 3 method [||. 
On comparing with experiment it is clear that there is a disagreement with theory on the sign of 
the spin density at the bridging carbon atoms. On should remember however, that the theory 
applies to an isolated cation in vacuum. Calculations for smaller than the formal charge reverse 
the spin density at the bridging carbons. The crystal environment may also change areas with 
small spin density. 



4.3 Spin susceptibility in the high-T c superconductor YBa2Cii307 x 

The discovery of the high-T c superconductor La 2 _ x Ba x Cu04 with T c = 35K by Bednorz 
and Miiller in 1986 has been at the origin of an enormous amount of work to under- 
stand the electronic (charge and spin) correlations in these materials. Following the discov- 
ery of the La 2 _ x Ba x Cu04 compound, other materials exhibiting similar or higher transition 
temperatures for superconductivity have been synthesized, like La 2 _ x Sr x Cu04,Nd 2 _ x Ce x Cu04, 
YBa 2 Cus07_ x , and others. All these materials share common features of their crystallographic 

1 1 



structure. They possess CuCVlayers well separated from each other, so that they can be consid- 
ered as quasi-two-dimensional materials. In this class of materials, superconductivity is achieved 
by carefully tuning the amount x of Sr, Ce or O which results in doping the CuCVlayers with 
charge carriers. The important feature is that the cuprate materials are either antiferromagnets 
and insulators or paramagnetic metals and superconductors below a critical temperature T c . 
Important antiferromagnetic correlations and fluctuations persist in the superconducting phase. 
The role played by these fluctuations in the formation of the superconducting state is still the 
subject of an intense debate. A central piece of the physics of the high-T c superconductors is 
the understanding of the charge and spin states in the CuCVlayers as a function of doping both 
below and above the transition temperature T c . 



The intensity of the scattered neutrons can be directly related to the spin susceptibility at Q = 
and uj = through the dissipation-fluctuation theorem. As the signal is particularly small in 
the high-T c compounds, use of polarisation analysis is required to enhance the contrast and to 
isolate the weak magnetic contribution. The results obtained for the temperature variation of 
the local susceptibility at the copper sites in the C11O2 layers is shown in Fig. |5|. A particularity 
of the temperature dependence of the signal is the appearance of the so-called spin pseudo-gap 
in under-doped YBa 2 Cu0 6 .52 for which the spin susceptibility drops above T c . On the other 
hand, the local spin susceptibility in optimal doped samples decreases only upon cooling below 
the superconducting temperature [M . 



Figure 6: Local spin susceptibility as a function of temperature on copper sites in the Cu0 2 
layers in the high-T c superconductors YBa2Cu307_ x (taken from Ref. PB[). 



5 Spherical Neutron Polarimetry 



Spherical Neutron Polarimetry (SNP) has recently been developed and successfully tested at 
the ILL ]50| as an alternative way of measuring magnetic structures. Moreover, this method 
allows to determine form factors and spin densities in antiferromagnets for which very few data 
is available. The classical technique discussed in chapter [| cannot be applied in antiferromag- 
nets with propagation vector Q = when the magnetic and nuclear structure factors are in 
phase quadrature. For such cases, like Cr 2 3 or even hematite, the neutron cross-section a is 
polarisation independent [[Hp 

P f a = P (l - 7 2 ) + 2 7 2 Q(Po ■ Q) + 2 7 (Po x Q), (17) 

with a = 1 + 7 2 and 7Q = 3D_i_(Q)/iV. SNP gives access to the complete set of independent 
correlation functions involved in the nuclear-magnetic scattering process by a direct measure- 
ment of the three components of the polarization vector P/ of the scattered neutrons. Eq. [5] 
shows that if the polarization Pj of the incoming neutron beam is fixed, a measurement of P/ 
allows in most cases an unambiguous determination of the direction of the magnetic interaction 
vector Dj_. This is an alternative way of determining magnetic structure factors to the standard 
diffraction method that relies on a precise measurement of neutron intensities. Measuring inten- 
sities is the same as measuring Dj_ • D^j_ which leads to a loss of phase factors and often magnetic 
structures cannot be unambiguously resolved by unpolarized neutron diffraction. SNP has been 
successfully applied in problems involving complex magnetic structures, like spiral structures, 
systems with magnetic domains and small magnetic moments, frustrated antiferromagnets, etc. 
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It has to be pointed out that the method is sensitive to the direction of the magnetic interaction 
vector only and not to its magnitude. In this case, finding the value of magnetic moments re- 
quires, as usual, the comparison of magnetic and nuclear cross-sections. In contrast to standard 
single-crystal diffraction, it is insensitive to secondary extinction and allows magnetic structure 
determination of samples even in the presence of magnetic domains. Namely, magnetic domains 
depolarize the neutron beam according to their respective population. In other words, the do- 
main population is obtained by measuring the amplitude of the neutron polarization vector 
Pf. In order to perform spherical neutron polarimetry a non-isotropic domain distribution is 
usually necessary. It can be induced for example by the application of uniaxial pressure. 



Spherical neutron polarimetry out-performs standard polarized neutron scattering as it allows 
to measure both the longitudinal and transverse components of Pf. Namely, if a magnetic 
field is applied to the sample, as it is the case for the longitudinal polarimetry explained in 
section ^, only the component of the polarization longitudinal to the field can be measured. The 



transverse components depolarize rapidly and are lost [52]. This is the case in antiferromagnets 



with mixed nuclear-magnetic Bragg reflections, where most information is contained in the 
transverse components through the nuclear-magnetic interference term. 



5.1 Realization of a zero-field chamber: Cryopad 

Following the introduction in the previous section it is clear that the transverse components of 
the polarization can only be measured if the sample is placed in a zero-field chamber. Such a 
device (called Cryopad) has been constructed at the ILL (Fig. [7]) []53| . It consists of two cylindric 
Meissner shields in the superconducting state. The diameter of the inner shield is large enough 
to accommodate a cryostat and/or other devices to define the sample environment. Cryopad 
is centered on the sample table and its orientation is fixed with respect to the wave-vector kj 
of the incident neutrons. Therefore, the sample can be oriented independently form Cryopad 
in order to access various Bragg-reflections. 

Figure 7: Schematic zero-field chamber Cryopad II used for spherical polarimetry at ILL (Taken 
fromRef. @). 

The components of the polarization of the incident and scattered neutrons are defined indepen- 
dently by means of two rotating solenoids (called nutators) that are placed in the incident and 
scattered neutron beam. They act as guide-fields in order to orient the neutron polarization 
vectors P a (a = incident, final) in the plane parallel to the Meissner shields (transverse to k a ). 
Two spin turning coils between the Meissner shields apply a horizontal field transverse to k Q 
and allow the definition of the component of P a along k Q . The modulus of P/ is determined 
by measuring the flipping ratio of the scattered neutrons simply by reversing the field of the 
nutator after the sample. This can be accomplished quickly by reversing the current in the 
solenoid. The combined use of the nutators and of the spin turning coils in Cryopad allows the 
analysis of P/ for any direction of P^. In practice, the information is obtained by measuring the 
three orthogonal components of the final polarization Pf (a = x, y, z) for the three orthogonal 
directions of the initial polarization Pj. The direction x is defined as being along the scattering 
vector Q, z is chosen perpendicular to the scattering plane and y is the last orthogonal direc- 
tion. For example, if the initial polarization is chosen along the z-axis, the component of the 
final polarization along the x-direction is given by P x = (n x+ — n x -)/(n x+ + n x -), where n x+ 
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and n x - are the number of neutrons with spins along or anti-parallel to the x-axis, respectively 
5"Bfl . Measuring the nine flipping ratios is sufficient to determine the required information about 



the orientation of the magnetic interaction vector Dj_(Q) and to obtain the value of the ratio 
between magnetic and nuclear amplitudes |Dj_(Q)|/iV. Geometrical relationships between the 
direction of the final neutron polarization relative to Pj have been derived by Nunez et al. []5i| 
and are of great help to determine the direction of D^(Q). 



5.2 Example: UPtGe 



Heavy-fermion materials are characterized by a very large linear coefficient of the specific heat 
and a greatly enhanced Pauli susceptibility, corresponding to effective masses of the quasi- 
particles of about two orders of magnitude larger than the free electron. Heavy fermions are 
therefore ideal systems to study strong electron correlations and the number of compounds 
showing heavy-fermion behavior is large. The ground-state of these systems varies from metal- 
lic (CeCu 6 , Celn 3 ) to insulating (Ce 3 Bi 4 Pt3) and from antiferromagnetic (U 2 Zn 17 , UPd 2 Al 3 , 
UNi 2 Al 3 ) to superconducting (UPt 3 , UPd 2 Al 3 , UNi 2 Al 3 , CeCu 2 Si 2 ). In most of these systems 
the magnetic ground state is determined by the competition between the Kondo-effect that tries 
to screen the magnetic moment and the RKKY interaction that tends to stabilize a ground-state 
with long-range magnetic order. 



Non-collinear magnetic structures in compounds with localized spin-densities can be explained 
on the basis of the Heisenberg model to originate from competition between exchange forces. 
In systems with 5/-electrons, like U 3 P4 or U 2 Pd 2 Sn it has been shown by calculations based 
on the local-spin density functional theory that a non-collinear arrangement of the magnetic 
moments is the consequence of strong spin-orbit coupling [p5 ]. However, this theory does not 
favor the helicoidal-type of magnetic structures found in e.g. UNi 2 Al 3 and UPtGe. 



As an example of a magnetic structure determination with the help of spherical neutron po- 
larimetry we show results obtained from the ternary compound UPtGe that orders below 
T N ~ 50 K with a propagation vector Q = (0.554,0,0) p6| . Measurements in single crys- 
tals using unpolarized neutrons could not decide between an amplitude-modulated spin-density 
wave and a cycloid with unequal magnetic moments along the a- and c-axis (ellipticity), respec- 
tively, yielding for the two models similar agreement factors between observed and calculated 
structure factors |57| . 



It is seen from Eq. |5] that the direction as well as the amplitude of the polarization of the 
scattered neutrons depend upon the value of the chiral term. In particular, the chiral component 
disappears, when D_|_(Q ) is parallel to D^(Q ), as it is the case for an amplitude modulated 
wave and the polarization of the neutron precesses by 180° around the scattering vector [Q. 
In contrast, the chiral term gives a contribution to P/ if the magnetic structure is a cycloid. 
Hence, the direction of P f depends on the scattering geometry. Spherical neutron polarimetry 
therefore allows to distinguish between an amplitude modulated spin-density wave and a helix 
H. For UPtGe, the direct ions of polarization of the diffracted neutrons for directions of Pj 
perpendicular to the scattering plane (z), along the scattering vector (x), and a third direction 
in the scattering plane but perpendicular to x are summarized in table |2|. The results of 
spherical neutron polarimetry unambiguously show that the magnetic structure of UPtGe is a 



cycloid with an axis ratio ~ 1.24 |58fl, as shown in Fig. 
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Figure 8: Magnetic structure of UPtGe as determined by spherical neutron polarimetry using 
the Cryopad device (after |58|1). 



Table 2: Spherical Neutron Polarimetry data obtained with Cryopad in UPtGe (after Ref. 
Pi, Pf and P ca ic are the incident, scattered and calculated neutron polarizations. The 
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polarization axis are defined as: x is parallel to the scattering vector Q; y is perpendicular to 
Q and in the scattering plane; z is vertical. 
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6 Inelastic Neutron Scattering with Polarized Neutrons 

In analogy to the neutron cross-section derived by Blume [fjj for the elastic case, there are three 
contributions to the inelastic cross section: 



• a pure nuclear one that gives rise to phonon scattering, 

• a pure magnetic one when neutrons are scattered e.g. by spin waves, 

• and a magnetic- nuclear interference term that is present only in special cases, as e.g. 
when the spin-lattice interaction in a ferromagnet cannot be neglected. 



The inelastic neutron cross-section and its relationship to the neutron polarisation have been 
derived by many authors (see e.g. ||T], |60]). In the following we will reproduce the calculations 
of Maleyev |6l| which expresses the time-dependent scattering amplitudes in terms of the Van 



Hove correlation function 

ttS ab (uj) = - — < A, B >», (18) 

1 — exp(— uo / 1 ) 

where < A, B >" is the absorptive part of the generalised retarded susceptibility 

< A,B >l= tt(1 - exp(-w/T))Z- 1 £ e x P (-^)A ab B ba 5(u + E ab ). (19) 

a,b 

A and B are operators like iV(Q) or Dj_(Q) from Eq. Q Z~ l is the partition function and 
E a ^ are the energies between eigenstates of the system. Expressed in such terms, the inelastic 



neutron cross-section is given by 

a = a n + a m + a nmi (20) 

with 

*n = 7 ^r<iV(-Q),iV(Q)>^ 

ki Ti 1 — exp(— LJ/T) 

* PI / J Q( g 7 

= 7 ^r(<iV(-Q),D ± (Q)>^ + <D ± (-Q),iV(Q)>:).P l . (21) 

ki til — exp{—Lu/l ) 

{a,/3,7} = {x, y, z} in Cartesian coordinates. The polarisation vector P/ is accordingly given 

• for nuclear scattering by P/cx n = Pi<7 n , 

• for pure magnetic scattering by 

PfaCTm = ^P l/3 (< ^(-Q),^(Q) >: + < ^(-Q),^(Q) >:) 

P 

- ]T <W < DK-Q), DKQ) >l -i J2 Co* < Di(-Q), DKQ) >» ut 

M Pi 

• for magnetic-nuclear interference scattering by 

P f a nm = <iV(-Q),D ± (Q)>: + <D ± (-Q),iV(Q)>^ 

+i{< iV(-Q),D ± (Q) >l-< D ± (-Q),JV(Q) >£] x P*. 

It turns out from these equations that whereas phonons do not change the neutron polarisation, 
scattering by spin waves do. As we will see below this feature is very useful to separate and 
identify the different magnetic modes in ferro- and antiferromagnets. For the particular cases 
of simple ferromagnets and of two-sublattices collinear antiferromagnets, explicit expressions 
for the polarisation dependence of the inelastic neutron cross-section can be found in the clas- 
sic paper of Izyumov and Maleev |62|. The magnetic-nuclear interference term gives rise in 



particular to the so-called magneto -vibrational scattering and is also important for a proper 
understanding of magneto- elastic scattering. Their origins are due to the fact that the cloud of 
electrons that carries magnetism follows the nuclei when they oscillate around their equilibrium 
position and that the magnetic moment is modulated by the lattice vibrations, respectively. 
The magneto- vibrational scattering is inelastic in the nuclear system but elastic in the magnetic 
one. It occurs at the same positions in reciprocal space as the phonons but with a polarisa- 
tion dependence. It has been exploited e.g. to measure the magnetic form factor through the 
polarisation dependence of phonons at general Q positions in Fe [p5[. 



6.1 Longitudinal neutron polarimetry 



The first spectrometer that allowed the analysis of the scattered neutrons was built by Moon, 
Riste and Koehler by replacing monochromator and analyzer of a triple-axis spectrometer by 



ferromagnetic crystals that were saturated in a magnetic field (Fig. |9|) Q . In contrast to spher- 
ical neutron polarimetry the polarization is maintained by means of guide fields throughout the 
instrument. As a consequence, only the component of the neutron polarization parallel or anti- 
parallel to the field direction can be measured whereas the transverse components depolarize 
and are lost. In exactly the same way than the magnetic-nuclear interference term cannot be 
measured by standard diffraction, the same effect happens for the inelastic counterpart. The 
magnetic-nuclear interference contribution leads to a rotation of the initial polarisation which 
averages out if a magnetic field is applied. Such a contribution is only accessible if the sample is 
enclosed in a zero-field chamber, like the Cryopad device. The standard triple-axis instrument 
with polarisation analysis developed by Moon et al. allows to measure the energy, momentum 
and spin dependence of cross sections with the restriction that only longitudinal polarimetry 
can be performed. In practice, two basic scattering geometries are commonly used: Namely 
the spin-flip and non spin-flip cross sections are measured with either the polarization of the 
neutrons parallel or perpendicular to the scattering vector Q. 



Figure 9: Schematic arrangement of the three axis spectrometer for polarized neutrons used by 
Moon et al. ||], at Oak Ridge National Laboratory. 

An immediate application of longitudinal polarimetry is that for P, || Q all magnetic scattering 
involves processes in which the spin of the neutrons is flipped because D± z = 0. We point out 
again (see Eq. |3p that this is a general rule valid for elastic and inelastic as well as coherent 
and incoherent scattering. In contrast, if the scattering geometry is chosen such that Pj _L Q 
and M || Pj then the elastic magnetic scattering (D x = D y = 0) is non spin-flip, whereas the 
inelastic scattering is spin-flip (5~D y , transverse excitations) and non spin-flip (<5D 2 , longitudinal 
excitations), respectively. 

Polarized neutrons in inelastic neutron scattering are often used to separate the magnetic from 
the nuclear scattering or to distinguish magnetic fluctuations perpendicular and transverse to 
the magnetization or scattering vector. When scattering by phonons dominates a neutron 
spectrum an unambiguous determination of the magnetic contribution to the neutron cross- 
section can be accomplished by measuring the scattering with P.; || Q. A typical example for 
UFe2 is shown in Fig. where the linear dispersion curve of the acoustic phonons can be 
distinguished from the quadratic dispersion curve of the spin-wave branch in UFe2 at low T in 
the cross-over regime [B4[ . 



Figure 10: Constant-energy inelastic scan in UFe2 using polarization analysis showing that 
spin waves occur in the spin-flip channel (black symbols), while scattering by phonons is non 
spin-flip (open symbols) (taken from Ref. |64[]). The line is simply to guide the eye. See text 
for details. 



Longitudinal polarimetry is not only an important method for measuring magnetic and nuclear 
cross sections unambiguously, it is also very powerful in separating self and collective dynamics 
in materials that contain strong incoherent scatterers like hydrogen in biological materials and 
polymers |)5| (see chapter [?]) . 



6.2 The XYZ method 



The longitudinal polarimetry method can be generalized to polarization analysis along the 
three Cartesian directions, the so called XYZ method (that is however to be distinguished from 
spherical polarimetry). This method allows to apply the technique of polarization analysis to 
time-of-flight spectrometers with multi-detectors. With the coordinate system shown in Fig. 



Ill the non-spin flip (NSF) and spin-flip (SF) cross- sect ions are given by |66| 
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Here, om is the magnetic, <tns is the nuclear spin incoherent and <tn is the coherent plus 
isotopic incoherent nuclear scattering cross-section. The magnetic scattering can be isolated 
by combining these equations in the following way (independent of the angle a) 



a M - 2a z NSF - a y NSF - a x NSF (28) 
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o* SF + a y SF - 2a z SF . (29) 



The XYZ-method has been applied successfully in determining the dynamical magnetic response 
in metals with small magnetic moments like V2O3 and in probing both atomic and spin 
correlations e.g. in spin-glasses @ or magnetic defects in disordered alloys. A full account of 



the applicability of the XYZ technique to this problem has been recently reviewed by Cywinski 
et al. |68|]. Another important application of the XYZ-method is that it allows to separate 
incoherent and coherent atomic motions, as presented in chapter 0. 



Figure 11: Geometry of the XYZ polarization method. Pj and P/ are the polarizations of the 
incident and scattered beams, respectively. Q is the scattering vector. 



6.3 Paramagnetic scattering 

According to the Rhodes and Wohlfarth J7D| theory, magnetic materials with ci-electrons can be 
classified in localized and itinerant systems. While for systems with localized spin densities, the 
magnetic moment in the paramagnetic phase is temperature independent, the ratio between 
paramagnetic and ordered moments varies with temperature in the Stoner model. The theory 
of spin fluctuations for localized and itinerant magnetic systems is reasonably well developed 
in the paramagnetic phase ||71|| . In that respect, inelastic scattering of neutrons provides direct 
experimental information on the spectrum of spin fluctuations on an absolute scale as it gives 
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access to the space- and time- variation of the spin-spin correlation function S(Q,u) that is 
related to the imaginary part of the dynamical susceptibility x(Q,cj) [|TJ. Paramagnetic scat- 
tering is usually very weak and difficult to separate from coherent (phonons) and incoherent 
nuclear scattering. However, the signal can be uniquely identified in experiments by using the 
difference method, namely, the difference between the spin-flip scattering as measured in a 



(small) field parallel to Q and perpendicular to Q contains only magnetic scattering |K| [72fl . 
The reason being that, the inelastic scattering from phonons is suppressed and nuclear inco- 
herent scattering as well as room background cancel. This statement is generally valid as long 
as the nuclear magnetic moments are disordered, i.e. (I 2 ) = {I 2 ) = (I 2 ) = |/(J + 1), and 
(I a ) = 0. Nuclear ordering occurs only for extremely low temperatures |73| . 



Therefore one obtains for example 
1 
2 

Once the intensity of the paramagnetic fluctuations is measured, the ^-integrated intensity can 
be put on an absolute scale by comparison with an acoustic phonon measured close to a Bragg 
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peak [|74| or by using an incoherent scatterer like vanadium. Hence, an effective, paramagnetic 
moment can be found, defined by 



M{q) 



/(?) 



S(q, ui)dui 



1/2 



(31) 



where f{q) is the form factor |72 



74j| pointed out that Eq. |3l] overestimates 



Ishikawa et al. 

the amplitude of spin fluctuations at low temperatures when the energy range of the spin 
fluctuations extends beyond ~ k^T c . They propose instead to use the Kramers-Kronig relation 
to obtain first the static susceptibility that is linked via the fluctuation-dissipation theorem to 
the amplitude of the spin fluctuations through (M 2 (q)) = 3fcsTx(g) for hu <C k B T and 

S(q, u)[l — exp(—hu/k B T)] 



2 2 
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-duo. 



(32) 



Measurements of the paramagnetic fluctuations in MnSi with a coarse energy resolution 
that the ^-integration is automatically performed) with Eq. |3l] |72] and Eq. |32 



so 

show that 



the amplitude of the local magnetic moment indeed increases with increasing temperature, in 



agreement with self-consistent renormalisation theory (Fig. |12D [[71 



Figure 12: Temperature dependence of 47rq 2 (Mq) plotted against q in MnSi as measured with 
polarized neutrons by Ishikawa et al. |[74|| . It is apparent from the figure that the mean-square 
amplitude of the spin fluctuations increases with increasing temperature in agreement with the 
calculations of the self-consistent renormalisation theory |7T| . 



Using the difference technique, the scaling behavior of many different itinerant ferromagnets 
has been investigated in the paramagnetic phase and it was shown that the scattering functions 
of Fe and Ni can be modeled above Tc by a simple Lorentzian scattering function given by J75] 

S(q, w) = — j7^x(q = o) ? K , 9 -n 9 rg 9 , (33) 

1 — exp {—uj/1 ) 
where the inverse correlation length k = k,q(T/T ( 

f(x) is approximately given by the Resibois-Piette scaling function [76 



h? + q 2 T 2 + uj 2 ' 
l) - 7 , the line-width T q = Aq 2 - 5 f(n/q), and 



To conclude, polarized neutron scattering is a powerful method to measure paramagnetic fluc- 
tuations in particular when the signal is weak. 
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6.4 Transverse and longitudinal excitations in ferromagnets 



The magnetic properties of compounds with localized spin densities are usually described by 
the Heisenberg Hamiltonian 

H=-Y^J ij S i -S j (34) 
hi 

where Jjj is the exchange integral between the spins located at the i- and j-position, respectively 



Tt\ . Depending on the sign of the exchange integral, Eq. |34| favors either antiferromagnetic 
or ferromagnetic ground-states. If exchange interactions extend beyond nearest neighbors, 
competing effects can occur that may lead to non-collinear or even incommensurate magnetic 
structures. 

Because of its simplicity, the Heisenberg ferromagnet is often taken as model system to study 
the properties of phase transitions. Within the simple picture of localized spins, long range 
order is lost due to the thermal excitation of spin waves that evolve into the critical scattering at 
To- Using unpolarized neutron scattering the spin dynamics close to Tq has been investigated 
in detail [fT8| |. It was shown that the spin waves, i.e. the spin fluctuations transverse to the 



magnetisation vector M, renormalise close to T c and that the susceptibility x{o) as measured 
at small angles diverges at small q for T — > Tq. Here, q is the reduced momentum transfer with 
respect to the nearest Bragg peak. The divergence of xil) is due to longitudinal fluctuations 
along M because the cross section for spin waves does not contain a correlation length that 
diverges at T C - Unpolarized neutron scattering was not successful in detecting the longitudinal 
fluctuations in ferromagnets in contrast to the situation in antiferromagnets J79L 80 1, where 
they can be easily observed. 

The longitudinal fluctuations can be isolated by means of inelastic neutron scattering with 
polarization analysis f81 |. The experiment is performed by measuring the differential spin-flip 



and non spin-flip cross sections from a ferromagnetic sample, for example EuS, that is saturated 



in a vertical magnetic field B„ that is perpendicular to the scattering vector Q. Fig. [L3| shows 
three typical measurements on EuS that have been performed longitudinal and transverse to 
the reciprocal lattice point (200) at 0.93Tc |jS2|. 



Figure 13: Constant-Q scans probing magnetic fluctuations in the ferromagnetic phase of EuS. 
The solid lines are fits to the data using Lorentzian spectral weight functions convoluted with 
the resolution function of the spectrometer IN14 at the ILL. The longitudinal spin waves are 
reduced in intensity due to the dipolar interactions. The parallel fluctuations are quasielastic. 

The spin-flip data shows the spin waves with polarization vectors 5S transverse and parallel 
to the reduced momentum transfer q. The former are the Goldstone modes of the system and 
diverge like x^ w oc 1/q 2 (Table Q). The longitudinal spin waves attain a mass |33 due to the 



dipolar interactions and do not diverge, Xsw ^ l/(<? 2 + <?d)- Finally the non spin-flip data 
shows the longitudinal fluctuations that are quasielastic and diverge like Xz oc l/(g 2 + k 2 ). 
Because the width Y q of Xz{o.i u ) ls comparable to the spin wave energy E q it is clear why the 
longitudinal fluctuations escaped detection with unpolarized neutron scattering. The results 
are in qualitative and quantitative agreement with a coupled mode analysis |S4 and mode-mode 
coupling theory MB. 
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Table 3: Transverse and longitudinal susceptibilities of a Heisenberg ferromagnet with dipolar 
interactions in the ordered and paramagnetic phases for different directions of the momentum- 
transfer q and accessible by one-dimensional polarization analysis (taken from Ref. ||82|| ). 
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g 2 q 2 + K 2 +q 2 D q 2 +K 2 q 2 q 2 +K 2 +q 2 
1,1,1 2,2, 1 
q 2 q 2 +q 2 n q 2 +K 2 q 2 +K 2 1 q 2 1 q 2 + K 2 +q 2 



6.5 Spin waves and phasons in incommensurate, antiferromagnetic 
Cr 



One of the outstanding features of antiferromagnetic Cr is the occurrence of an incommensurate 
spin-density wave below T N = 311 K that is transversely polarized (S perpendicular Q ± ) 
with Q ± = (1 ± 5,0,0) being the incommensurate wave- vector f86f . The magnetic excitations 
exhibit many unusual features that are not well understood. In particular, the magnetic modes 
that originate from the magnetic satellite peaks at Q 1 * 1 have such a steep dispersion that the 
creation and annihilation peaks cannot be resolved anymore. Using inelastic scattering of 
unpolarized neutrons and analyzing the width of the peaks in constant energy scans led to the 
conclusion that the velocity of the excitations is c sw ~ 1020 meVA @. This value deviates 



significantly from the theoretical value of a random phase approximation (RPA) that is given 



by = y ~vf — 1500 meVA, where vp is the Fermi velocity |||. Using unpolarized neutron 

scattering it has been shown that transverse (with respect to the staggered magnetization) as 
well as longitudinal excitations contribute to the inelastic scattering that emerges from the 
incommensurate satellite peaks pi . 



In the absence of sizable magnetic-nuclear interference contributions in the cross-section, the 
most direct way to separate the transverse from the longitudinal fluctuations is the use of 
longitudinal polarimetry. For such an experiment it is necessary to use a Cr single-crystal 
cooled through T N in a large magnetic field in order to induce a single-Q state. During the 
experiment, a vertical field B = 4 T was applied along [001] in order to enforce a single- 
domain spin density wave with the magnetic moments aligned along the [010] direction. In this 
configuration, the spin-flip scattering is due to the longitudinal modes and the non spin-flip 
scattering due to the transverse modes. 



Fig. [14] shows constant energy scans for E = 4.2 meV measured in the transverse spin-density- 
wave phase at T = 230 K (0.74T/v) |JIJ. It is clearly seen that the inelastic, incommensurate 
peaks with transverse polarization are significantly sharper than the corresponding longitudinal 
peaks. Therefore, the mode velocity of the spin waves, c sw , is significantly larger than the mode 
velocity of the phason modes, c p h- These results are in qualitative agreement with results of 



RPA theory [pqj . In addition, the data shows that the enhanced magnetic scattering at (100) 
and E = 4.2 meV has a longitudinal polarization. Without going into further details the results 
indicate that a proper understanding of the magnetic excitations in Cr can only be gained if 
polarization analysis is used. 
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Figure 14: Constant-i? scans at 4.2 meV, probing the longitudinal and transverse excitations 
along the [100] direction in the transverse spin- density-wave phase of Cr at T = 230 K. The 
inset shows the intersection of the constant-i? scan with the dispersion of the the transverse 
(solid lines) and the longitudinal modes (broken lines) as well as the Fincher-Burke modes. 

6.6 Magnetic excitations in a heavy fermion superconductor 



The heavy fermion superconductor UPd2Al 3 exhibits the unusual coexistence of antiferromag- 
netism and superconductivity below T c = 2 K, i.e. the ordered magnetic moments of the f- 
electrons of U persist in the superconducting phase [51] . This has been taken as a sign that the 
interplay of magnetism and superconductivity could be studied in this material. Neutron p2 
and x-ray scattering |93] experiments have shown that the magnetic structure of UPd2Ai3 con- 
sists of ferromagnetic planes stacked along the c-axis with a propagation vector Q = (0, 0, 0.5). 
The magnetic moments are confined within the hexagonal plane and are found to have an un- 
usually large value of /i = 0.85/iB at saturation. First elastic |95[ and inelastic |96| neutron 
scattering experiments could not unambiguously reveal any change in the magnetic properties 
of UPd2Al 3 upon cooling the sample below the superconducting transition temperature. 



Inelastic neutron scattering experiments |97|, |98| performed with an improved energy-resolution 
as compared to the work of Petersen et al. |96|] showed that there exist two contributions 
to the spectrum of magnetic fluctuations in UPd2Al3. While the first one corresponds to the 
spin- wave previously measured by Petersen et al. [ |56[] , a second mode localized around the 
antiferromagnetic wave- vector Q is observed in the energy range < E < 0.5 meV. This lower 
energy mode is heavily damped for all temperatures in the antiferromagnetically ordered phase 
and strongly sharpens upon passing into the superconducting phase. At the lowest temperature 
the low-energy mode develops an apparent energy gap with a value comparable to T c |)8|, f39| . 
Further evidence of a strong interplay between magnetic fluctuations and superconductivity 
in this compound originates from the use of polarized neutrons as shown in Fig. [y| Using 
a polarized beam, it was possible to show that the two magnetic modes are both polarized 
transverse to the magnetization vector and hence are likely to interact with each other ||100|| . 

To perform this experiment the sample was field-cooled, so that the magnetic domains could 
be aligned along the magnetization vector M. Using a neutron polarization perpendicular to 
the scattering plane, it turns out that magnetic fluctuations parallel to the magnetization are 
non-spin flip, while those perpendicular to M appear in the spin-flip channel. Analysis of the 
line-shape of the inelastic neutron scattering data suggests that /-electrons located in a small 
energy range around the Fermi surface play a significant role in forming the superconducting 
state in UPd 2 Al 3 fllOl 



Figure 15: Experimental data from UPd2Al 3 at the antiferromagnetic wave vector Q = 
(0, 0, 0.5) and T = 150 mK. The data were taken with a fixed outgoing neutron wave vec- 
tor of kf = 1.15 A -1 . Frame a) and b) are taken with polarized neutrons. For means of 
comparison a scan measured with unpolarized neutrons is shown in frame c). In frame a) and 
b), the transverse response is shown as black symbols, whereas the longitudinal component is 
represented by open circles. See Ref. ||100|1 for details. 



6.7 Magnons and Solitons in low-dimensional systems 



The magnetic properties of low- dimensional compounds have attracted a lot of attention as 
new effects due to quantum fluctuations are strong. For one-dimensional Heisenberg antifer- 
romagnets the ground states and energy excitations are different for integer- and half-integer 
spins ||102|| . Antiferromagnetic chains with S —1/2 spins have a disordered ground-state. The 
low-lying excitations are characterized by a continuum of excitations without energy gap at the 
zone center. On the other hand, for integer-spins a finite energy gap was predicted by Haldane 
and obtained by numerical calculations ||102| , |103|| . Examples of materials exhibiting a Haldane 



gap are NENP [|04 



The characteristics expected for a 



Y 2 BaNi0 5 HI, or CsNiCl 3 
Haldane system have been observed in these compounds by inelastic neutron scattering, like i) 
a periodicity of 2tt in the magnon dispersion, ii) line-width broadening of the magnetic excita- 
tions as a function of momentum transfer indicating the presence of a two-magnon continuum, 
and iii) a large field dependence of the magnetic excitations 
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Polarized neutron scattering experiments have shown that the energy gap in the spectrum of 
magnetic excitations in CsNiCFj is a triplet ||108 . In such quasi-one-dimensional antiferromag- 



nets, weak inter-chain exchange interactions J' can lead to a Neel phase at low temperatures. 
In fact the ordering temperature depends on the ratio of the intra-chain interactions J to J' . 
Interestingly, for such systems, where antiferromagnetic ordering is close to disorder, linear 
spin-wave theory does not account properly neither for the energy-dependence of the magnetic 
excitations nor for the number of magnetic modes. In particular, the existence of a longitu- 
dinally polarized magnetic mode that cannot be predicted by standard or modified spin-wave 
theory has been proven by means of inelastic polarized neutron scattering in CsNiCl3 [|106|| and 
Nd2BaNiC>5 (Fig. |TBp | |109| ]. The instrumental set-up was chosen so that magnetic excitations 



transverse to the magnetic moments could be separated from the fluctuations along the spin 
direction. With such a geometry, it could be shown that additional excitations with longitudi- 
nal polarization are present in the spectrum of S(Q,u), in agreement with calculations based 
on renormalization-group theory jlT0|. 



Figure 16: Temperature dependence of constant Q-scans measured in Nd2BaNiOs with polar- 
ized neutrons. Open and black circles refer to spin-flip and non spin-flip scattering respectively 
(Taken from Ref. [TTUg]). 

The quasi-one-dimensional 5 = 1/2 inorganic compound CuGeOs presents the particularity to 
undergo a chemical phase transition below T = 14 K to a phase, called Spin-Peierls phase, 
where the copper chain is dimerised ||111|| . Consequently, the exchange interactions along the 



chain direction are not uniform anymore but alternate with values J and J', respectively. The 
magneto-elastic interaction is presumably responsible for this phase transition characterized by 
a non-magnetic ground state. For such a system, the spectrum of magnetic excitations attains 
a gap at the zone center whereas the first excited states are triplets. In CuGeOa the gap has 
a value of A ~ 2.5 meV ||112|| while away from the zone center, the spectrum of magnetic 
fluctuations is strongly dispersive along the copper chain direction. High-resolution inelastic 
polarized-neutron experiments, however, revealed that there is a second energy gap in this 
compound which separates the low-energy magnon-like mode from a continuum of excitations 
extending to higher energies ||113|| . The occurrence of two energy gaps in the spectrum of 



magnetic excitations in CuGeOs is clearly a signature of strong quantum fluctuations in S = 1/2 
antiferromagnetic chains. This behavior differs drastically from one- dimensional systems with 



large spin number S that in some cases can be described by the classical sine-Gordon equation. 
The combined effects of non-linearity and dispersion in these systems lead in addition to the 
linear excitations to a special class of excitations called 'solitons' ||108| . 



In a ferromagnetic chain, an excitation of soliton-type can be viewed as a 2ir turn of the spins 
over a small distance, in a similar way to a domain wall that would propagate through the 
crystal. These excitations are accessible to inelastic neutron scattering and in particular to 
polarized neutrons that allow to measure selectively the different space and time correlation 
functions S X ' X (Q, u), S y ' y (Q,u) and S z ' z (Q,u) ||108|| . Although experiments with unpolarized 



neutrons have shown strong evidence of non-linear excitations in the chain compounds CsNiF 3 
and TMMC, a crucial test for the existence of solitonic excitations is the measurement of the 



components fluctuating along and perpendicular to an applied magnetic field [ |114|| 



By separating the longitudinal part ^"(Q, tu) of the dynamic structure factor from the trans 
verse part S'- l (Q,c<j) with longitudinal polarimetry analysis (see Fig. [FT]), Boucher et al. ||115 



were able to study the wave-vector and energy dependence of solitonic fluctuations in the an- 
tiferromagnetic chain compound (CDs^NMnCl^TMMC). This study lead to the result that 
amplitude and lifetime of the solitons are strongly affected by collisions with magnons and by 
mutual interactions. As a consequence, the line shape and the line width of the experimen- 
tal dynamical susceptibility differ from the actual theoretical calculations based upon the low 



density non-interacting soliton gas model [|108 . 



Figure 17: Spectra of magnetic excitations measured in TMMC with inelastic polarized neutron 
scattering showing that the transverse and longitudinal fluctuations are different. See text and 



Ref. ||115|| for details. 



7 Self and collective diffusive atomic motions 



Collective motions of light atoms in metals consist of two different processes [ |116|| . The first 
one can be viewed as pure diffusion of ions through the lattice while the second process involves 
cooperative hopping of mutually interacting particles. Hence, the neutron scattering functions 
contain incoherent and coherent scattering contributions that are given within the random 
phase approximation by Lorentzian functions centered around zero energy transfer 

q (n \ 1 D & 2 



tt(AQ 2 ) 2 + u; 2 
S(Q) D C Q 2 

Scok{Q,<*)= ^ {DcQ 2y +UJ 2- ( 35 ) 

S(Q) is the static structure factor and D t and D c are the coefficients of incoherent and coherent 
diffusion, respectively. 

To separate the two quasi-elastic scattering processes which appear simultaneously in the neu- 
tron spectrum, it is best to use polarization analysis. In analogy to paramagnetic scattering, 
coherent and incoherent processes can be isolated by calculating the difference between non 
spin-flip and spin-flip scattering. In the case of different isotopes and disordered nuclear spins, 



the matrix elements for non spin-flip and spin-flip scattering are given for the coherent cross- 
section by 



Coherent = <T = (N)^ (36) 
Coherent = °~ + = ( 37 ) 



and for the incoherent scattering by pO 



incoherent = a incoherent = (N 2 )iso ~ (N) \ s0 + - (B 2 I(I + l)) iso (38) 

2 2 

° incoherent = ° incoherent = g + "0 ) ™o-> (39) 



where (■ ■ -)i SO refers to isotopic averaging and / to the nuclear spin [ 60fl . In compounds that 
contain scatterers with one isotope only, the coherent cross-section is obtained by dividing the 
spin-flip scattering by 2 and subtracting the result from the non-spin flip intensity. In this case, 



all the incoherent scattering is spin-flip scattering. Fig. ITS] shows the results of measurements 



of the self and collective dynamics of deuterium in a single crystal of Nb using the time-of- 
flight spectrometer D7 at ILL with the method described in chapter p.2| p5|. It is clear that 



polarization analysis gives an unambiguous separation of the incoherent and coherent quasi- 
elastic signals over a large range of momentum transfers Q. Such a study allows to determine 
distance and direction of jump processes, hopping and residential times through the analysis 
of the Q-dependence and energy width of the Lorentzian functions in Eqs. 

Figure 18: Coherent and incoherent scattering processes observed and calculated in a'-NbD .7 
by means of polarized neutron scattering on the multi-counter time-of-flight spectrometer D7 
at the ILL (taken from Ref. f6~5|). 



8 Conclusions 



The examples given in the previous chapters have shown that neutron scattering with polarized 
neutrons has become a very important means to measure magnetic properties over a wide 
range of Q and u and to distinguish between coherent and incoherent excitations in materials. 
Most experiments with polarized neutrons are being performed up to now at sources providing 
continuous neutron beams. The main reason being that most neutron polarizers are more 
ideally suited for applications with constant wavelength. The recent advances in the field of 
supermirrors and 3 He filters have improved the situation. The new devices allow to extend 
polarization analysis to high neutron energies and to use large area detectors. Therefore, we 
expect that polarization analysis will also soon become a standard technique at pulsed neutron 
sources. 

Recently, new developments in the field of magnetism have emerged that rely strongly on 
new developments in polarized neutron scattering. As a first example, we mention systems 
that can be characterized by sets of exponents, that differ according to the space- and spin- 
dimensionality and hence can be grouped into universality classes. In frustrated spin systems, 
the order parameter includes a term describing the spin chirality C = [Si x S2]. A direct 
observation of the fluctuations of the chiral variable is, however, impossible with unpolarized 
neutron scattering as these are related to four-spin correlation functions. Because the chiral part 
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of the neutron cross-section is polarization dependent ||117|| it can be observed with polarization 
analysis. 



As a second example, we mention that the Cryopad technique opens the possibility to apply 
spherical neutron polarimetry in inelastic neutron scattering as it allows to measure the nuclear- 
magnetic interference term directly. Recently, Maleyev |JT| has reconsidered the implication 
of the nuclear-magnetic interference term (NMIT) for inelastic scattering and has shown that 
in analogy with elastic scattering, it leads to a dependence of the neutron cross-section upon 
Pj, namely to a finite polarization of the scattered neutrons and to a rotation of the initial 
polarization. In particular, Maleyev has shown that the part of the dynamical susceptibility 
due to the NMIT is non-zero if there is a spin-lattice interaction characterized by an axial vector, 
as it is the case e.g. for the Dzialoshinskii-Moriya (DM) interaction D • [Sj x Sj]. Indications of 
the importance of the DM-interaction in the spin-lattice coupling and hence on its accessibility 
through inelastic spherical neutron polarimetry originates from recent inelastic experiments 
performed in the spin-Peierls compound CuGeC>3 where a rotation of the final polarization P f 
has been detected when the incident polarization is chosen parallel to the scattering vector 



118|| . Such measurements have shown that spherical neutron polarimetry can be applied in the 
field of inelastic neutron scattering, although counting times are long due to poor statistics. A 
further application may be the study of magneto-elastic coupling that plays an important role 



in invar alloys ||119|| . Together with new theoretical interest on such problems, it is probable 



that this method will contribute to an improved understanding of phase transitions, where spin- 
lattice interactions are important, as such effects cannot be studied by standard polarization 
analysis. 
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